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Abstract. We estimate the effects of low temporal frequency modes in the time stream on sky maps such as 
expected from the PLANCK experiment - a satellite mission designed to image the sky in the microwave band. We 
perform the computations in a semi-analytic way based on a simple model of PLANCK observations, which permits 
an insight into the structure of noise correlations of PLANCK-like maps, without doing exact, computationally 
intensive numerical calculations. We show that, for a set of plausible scanning strategies, marginalization over 
temporal frequency modes with frequencies lower than the spin frequency of the satellite (~ 1/60 Hz) causes a 
nearly negligible deterioration of a quality of the resulting sky maps. We point out that this observation implies 
that it should be possible to successfully remove effects of long-term time domain parasitic signals from the 
PLANCK maps during the data analysis stage. 
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1. Introduction and motivation. 

The planck satellite 1 is designed to measure the mi- 
crowave sky with unprecedented sensitivity and angular 
resolution. While its primary objective is to character- 
ize the anisotropics in the Cosmic Microwave Background 
(CMB), the full sky maps in each of the 9 frequency chan- 
nels will be the true PLANCK legacy. 

However, making maps from an instrument like 
PLANCK is non-trivial, due to the way the sky is scanned 
and the long-term correlations in the noise properties 
of the instruments (so-called 1// noise). Though a lot 
of attention has been given to the long term stability 
of the PLANCK instrument during its design phase it is 
practically unavoidable that slow secular drifts and semi- 
periodic signals will be present in the actual planck data 
(e.g., Seiffert et al, 2001, Mennella et al., 2002) and will 
require a software solution. In fact long term drifts, on 
time scales of hours and longer, are conspicuous in the re- 
cent data of the WMAP satellite and have been found to 
be of sufficient importance that a special treatment was 
applied to the time ordered data prior to turning them 
into the sky maps (Hinshaw et al., 2003). It seems there- 
fore prudent to assume that some kind of an analogous 
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procedure will be required to minimize the impact of such 
effects on the planck maps. 

Here we attempt to estimate the magnitude of this im- 
pact and elucidate the interplay of those parasitic signals 
with the usual 1// low frequency tail of the time domain 
noise as expected for PLANCK. More specifically we inves- 
tigate two (related) questions. First we strive to gain an 
understanding of how much the 'striping' in a map is con- 
strained by low-frequency information in the time stream 
(and therefore liable to be compromised by the parasitic 
signal mentioned above) and how much it is constrained 
by intersecting scan patterns on the sky. We shall find 
that the presence of 1// noise at any realistic level means 
that almost all of our ability to control the low frequency 
modes comes from the overlap of scan paths in the spatial 
domain, rather than from long time-scale information in 
the time stream. Secondly we wish to gain an understand- 
ing of how parasitic signals or specific removal techniques 
affect the structure of the map error matrix. Though for 
definiteness we will focus hereafter on very specific prob- 
lems, the analysis presented should be useful in guiding 
considerations of other similar issues for planck as well 
as a variety of different experimental setups. 

The structure of this paper is as follows. We start with 
a quick review of map making techniques (Section 2) and 
present briefly our assumptions about planck, its scan- 
ning strategy and anticipated performance (Section 3). In 
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Section 4, we analyze the low temporal frequency problem 
adopting a time domain perspective, while in Section 5 we 
derive corresponding pixel domain constraints. These two 
sections are quite technical, and the reader may wish to 
skim these on a first pass. We compare both results in 
Section 6 and finally draw some conclusions in Section 7. 

2. Map making review 

A reconstruction of a map of the microwave sky from a 
sequence of measured time ordered data samples is by 
now a well understood problem (e.g., Wright, Hinshaw, 
& Bennett, 1996, Tegmark, 1997a, Stompor et al., 2002). 
An entire suite of approaches, ranging from a simple bin- 
ning of observations into pixels on the sky, through de- 
striping techniques to optimal minimum variance map- 
making, has been to date successfully implemented, in- 
vestigated in detail (e.g., Wright et al., 1996, Tegmark, 
1997b, Delabrouille, 1998, Borrill et al., 2000, Natoli et 
al., 2001, Dore et al., 2001, Maino et al., 2002) and in a 
handful of cases converted into publicly available software 
packages (e.g., MADCAP 2 , MAPCUMBA 3 , MADmap 4 ). 
These different approaches trade, to various degrees, sim- 
plicity and numerical speed for accuracy of the produced 
maps, and usually they can be derived either as a special 
case of, or an approximation to, a more general formalism 
based on a maximum likelihood approach to map-making. 
The latter provides a handy and concise algebraic formu- 
lation of the entire problem (Tegmark, 1996). It has been 
demonstrated within such a formalism that a number of 
relatively straightforward generalizations are possible al- 
lowing us, in a statistically strict manner, to account for a 
variety of systematic problems commonly troubling actual 
experimental data (Wright et al., 1996, Oliveira-daCosta 
et al., 1999, Stompor et al., 2002, Hinshaw et al., 2003). 
All these developments highlight the practical feasibility, 
flexibility and robustness of this approach. In this con- 
text a major task in a successful accomplishment of the 
map-making stage for any CMB experiment becomes the 
detection and characterization of systematic contributions 
plaguing real data sets. The problem is made more diffi- 
cult by the sheer size and complexity of current and future 
data sets. Consequently, for map-making practitioners, the 
map-making issue is far from being definitely resolved and 
keeps providing challenges on a case- by-case basis. 

In a standard map-making procedure (Tegmark 1997a, 
Borrill et al., 2000) a final uncertainty of the map is en- 
coded in a pixel-pixel noise correlation matrix. In the case 
of a Gaussian instrumental noise, this matrix provides a 
complete description of the map error. However, in most of 
circumstances of the interest, given the size of forthcom- 
ing maps, the computation and storing of such a matrix 
is prohibitive, let alone any investigation of its structure 
and its dependence on a presence of parasitic signals in 

2 http://crd.lbl.gov/~borrill 

3 http : //ulysse . iap. f r/ cmbsof t/mapcumba/ 
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the data, or on details of a removal technique applied to 
such a systematic. In such cases an understanding of a role 
of a particular systematic effect and its impact on the fi- 
nal map needs to be investigated by different usually more 
indirect and simplified means. 

3. planck - basic features. 

To elucidate the impact of instrumental parasitics on low 
temporal frequencies (generally understood hereafter as 
those lower than the satellite spin frequency) we construct 
a simple model of planck which captures all of the essen- 
tial features for our purposes. We assume that the satel- 
lite spins on its axis roughly once per minute (i S pin = 
l//spin = lmin), with any given detector beam sweep- 
ing out a circle (or more generally - following Wandelt 
& Gorski, 2001, and Wandelt & Hansen, 2003 - a basic 
scan path) on the sky with opening angle of 29 = 170°. 
We will refer to a short stretch of a circle (ring) of a 
beam-size length as a ring pixel. The sampling rate is as- 
sumed to be A s = (^fwhm/3)/(11800/ sp i n ), corresponding 
to three measurements per the FWHM of a detector beam 
(= #fwhm) at each passing. Hereafter, whenever needed we 
will assume that #f w hm = 10'. Each circle is observed for 
lhour, before the satellite axis re-points and a different 
ring (shifted with respect to the previous one by ~ 2.5' 
along the great circle on the sky) is observed for another 
time T. We will assume that the re-pointing from a cir- 
cle to a next one is always instantaneous. As the satellite 
axis is re-pointed only by ~ 2.5' every hour, during every 
four hour long period each detector of the satellite ob- 
serves approximately the same sky sweeping during that 
time a single ring of a beam-size width. Hence, hereafter 
as basic scan paths we will consider sky rings as swept 
out during a time period of T = 4hours. The re-pointing 
frequency is thus, / rpt = 1/T ~ 7 • 10~ 5 Hz. The number 
of scans of each ring is T/i sp ; n and there are n rp i x (beam 
size) ring pixels per ring. As we will discuss in some detail 
in Sect. 5.1, n rp i x also corresponds to a number of (nearly) 
uncorrelated, independent ring pixels for a single ring of 
a PLANCK-like experiment and is therefore an important 
factor in the scaling relations of the pixel domain noise 
properties derived in the following. 

We should emphasize that though the discussion pre- 
sented here touches upon and depends upon a number of 
scan-related assumptions, the issues involved in optimiz- 
ing the scanning strategy for PLANCK are more complex 
and varied than what is included in the discussion pre- 
sented below. Clearly a more pointed analysis is required 
to elaborate on those issues. We leave such a discussion 
for future work. As far as this paper is concerned, while it 
is clear the numerical prefactors in our results are depen- 
dent on a particular scanning strategy, we believe that our 
main conclusions are independent on their specific details. 

Hereafter, we assume that the Gaussian instrumental 
noise can be accurately described as, 

P(/t) = <7 t 2 A s + (1) 
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where ft is a temporal frequency, /knee and at 2 parameter- 
ize the spectrum, and, a t describes the (white) noise level 
per measurement. Both these values are strongly detector- 
dependent. In the following, whenever necessary, we will 
assume the following values: for hfi, /knee — 30mHz 
and a t ~ lOOO^K (~ WOfxK^/sec/ y/S^) and for LFI, 
/ knee ~ 50mHz and er t ~ 3000,uK (~ 300^KVsec/ v / A s ~) 
(in both cases in antenna units). Here, hfi and lfi refer 
to the High and Low Frequency Instruments of PLANCK 
respectively. The actual flight parameters may be some- 
what better. However, for any anticipated values of /k nee 
the following relation is satisfied, 

/rpt ^ /spin ^ /knee- (2) 

As it is discussed below, our conclusions will not be af- 
fected by changes of either the noise level, which will only 
lead to a rescaling of the estimated values, or the assumed 
values of /knee as long as these are larger than £ lmHz. 

4. Time domain considerations. 

4.1. Background 

The generalized least squares map-making starts by mod- 
eling the observation as 

t = As + n (3) 

where t indicates the time stream data, A is the pointing 
matrix, s the sky signal and n the noise. If At = (ran*) is 
the (time-domain) noise correlation matrix, the minimum 
variance estimator of the sky signal, to, is 

m = [A*Aff 1 A] " 1 AW t ~ H (4) 

with [ J 4*A r t T 1 A] 1 the pixel-pixel noise correlation ma- 
trix. This equation can be efficiently solved using itera- 
tive methods if At is known (Dore et al., 2001, Natoli et 
al., 2001). Those methods do not require explicit com- 
putations of the pixel-pixel noise correlation matrix - a 
property which makes solving for the map estimate feasi- 
ble from both CPU time and disk/memory storage points 
of view. 

Clearly, in this equation two kinds of constraints are 
incorporated: first based on the known correlations in the 
time domain, as given by the noise correlation matrix, At , 
and second based on the fact that some parts of the (un- 
changeable by definition) sky are observed multiple times 
during an experiment, as encoded in the pointing matrix, 
A. Both constraints are usually tightly combined together, 
e.g., as seen in Eq. 4, giving us the 'best' possible map es- 
timate, for instance, that with a minimum noise variance. 

Efforts to remove/minimize systematic effects present 
on the time stream level can affect both of these types 
of constraints to a different extent. If, for instance, entire 
stretches of the time ordered data stream are rendered un- 
usable and effectively removed from the data, some pixels 
on the sky may not be observed multiple times anymore or 
even observed at all, clearly affecting the strength of the 



pixel domain constraints. If, in the contrary, only some 
of the temporal frequency bands have been found to be 
compromised by a systematic effect, it will be time do- 
main constraints which will be affected more directly. 

In the following we will focus on the latter case assum- 
ing that only low temporal frequencies have been contam- 
inated by a systematic effect of some kind but no time 
samples have had to be removed and therefore no changes 
to an actual scan pattern observed on the sky have been 
made. 

4.2. Destriping vs optimal map-making for planck. 

Given the rather simple scanning strategies planned for 
PLANCK, the sky map recovery can proceed in a number 
of ways. In particular, multiple re-scanning of the succes- 
sive rings on the sky opens up the possibility of performing 
map-making in two stages. In this case, first maps of the 
single rings are created which are subsequently combined 
together using the fact that different rings cross each other 
multiple times on the sky, i.e., exploiting the pixel domain 
constraint in a parlance of the previous Section. Such a 
two-step approach is a keystone of the destriping methods 
(Dcllabrouile, 1998, Maino et al., 2002, Keihanen et al., 
2003). Techniques of this kind are in general not optimal. 
That is the case even if single ring maps are produced us- 
ing optimal map-making (as we will assume in this paper 
henceforth), rather than simple sky binning as in 'tradi- 
tional' destriping. This fact can be understood as follows. 
As long as the beam is simply sweeping the sky in a cir- 
cle, in the frequency domain the sky signal is concentrated 
only at harmonics of the spin period and at zero frequency. 
The latter contains not only actual sky monopole but also 
all the higher multipole moments which do not average 
to zero over the circle. All together, these contributions 
amount to a constant ring offset. In a total power CMB 
experiment the constant offset of a single segment of the 
time ordered data is usually lost, leading to a loss of the 
ring offsets together with the corresponding part of the 
actual sky signal. 

Let us assume now that after spending time T on a 
given ring, the observation continues and the instrument 
is re-pointed and a scan of another circle commences, as 
it is the case for PLANCK. In the destriping methods, the 
continuity of the time stream from a ring to ring is ig- 
nored, and each ring is analyzed separately. Consequently, 
as in a single ring case mentioned above, the offsets of all 
the rings are lost during the initial single-ring processing 
and need to be subsequently restored (from now on only 
relative). The uncertainty involved in such a procedure in- 
creases the effective noise level in the map and a noisier 
map results. The situation looks different in the methods 
in which the time stream is not cut into segments, which 
are then analyzed separately, but treat the time stream in 
its entirety. Indeed, if the re-pointing is performed with- 
out interrupting the time stream continuity, the sky signal 
formerly confined to the zero frequency mode and the spin 
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frequency harmonics, now partially resides at the harmon- 
ics of the re-pointing frequency / rp t(= and at side- 
bands of the harmonics of the spin frequency (i.e., / = 

/spin ^ /rpt , /spin ^ 2 /rpt ; ■ ■ ■ , 2/gpmi / r pt , 2/ S pi n i 2/ r pt , . . . , 

etc) as well. Though the zero frequency mode is again lost, 
all the signal confined to the non-zero frequency bands is 
in principle accessible, rendering tighter constraints on the 
recovered sky map. If the instrument is successively and 
smoothly re-pointed in such a way that the observed sky 
area progressively increases, more of the multipole mo- 
ments of the map are contained in non-zero frequencies. 
In the case of a full sky covered in such a way, only the sky 
monopolc would be lost together with the zero frequency 
mode. 

In the following we address the issue of how much of 
an actual improvement over the two step method can be 
expected from such an approach in a case of a realistic 
PLANCK-like experiment. Given that, on the ring level, 
the two-stage method exploits only the pixel domain con- 
straints, while the full optimal map-making attempts to 
make use of both pixel and time domain ones, this question 
is clearly just a rephrasing of the problem posed in the pre- 
vious Section. At the computational level we can turn this 
problem into a question of how the time domain derived 
constraints on the relative offsets of two rings compare 
with those derived from the pixel domain. Alternately, as 
uncertainty in the recovery of the rings offsets results in 
the presence of strongly correlated linear features in the 
map aligned with the scan direction (i.e., stripes), we will 
refer to that problem as striping. 

The attractive advantage of the two-step map-making 
is its robustness to low frequency parasitic signal removal. 
Indeed, if the parasitic signal is confined to the frequency 
range / < / mar g ~ /spin, then high-pass filtering that part 
of the time stream which comes from each ring can remove 
most non-cosmological signals and at most an offset in the 
sky signal, causing no extra loss in precision. By contrast, 
the loss of the precision of the optimal map-making in- 
curred due to a removal of the contaminated low frequency 
modes needs to be investigated in detail (see Section 4.3). 

At this point it is worth emphasizing two issues. First, 
in the case of PLANCK, where /knee > /spin, if the final so- 
lution is to be nearly optimal it may need to account for 
the presence of I // noise within each ring. This suggests 
that Eq. 4 be solved on each ring (see Sect. 5.1) as the 
first step of the destriping method. Second, in this paper 
we will use the words filtering and marginalization over 
a given frequency band interchangeably, always meaning 
the latter. Henceforth, marginalization (or filtering) is un- 
derstood as a procedure of "weighting-out" unwanted fre- 
quency modes from the final map and the corresponding 
map noise correlation matrix by setting the noise level in 
those modes to infinity (i.e., the weights to zero). In both 
cases the relevant operations are usually done in Fourier 
space. More details on this approach can be found in the 
next Section and elsewhere (e.g., Stompor et al., 2002). 

Of course, the marginalization over the frequency band 
is not the only way of dealing with the parasitic signal. 



Other options can be viable, especially if extra information 
about the origin or character of the problem is available 
(e.g., Stompor et al. 2002). By comparison the marginal- 
ization may look like quite a drastic approach. However, 
as we show that in the following, for the low temporal fre- 
quency parasitic signal the marginalization turns out to 
be as good an approach as any other and more generally 
applicable. 

4.3. Constraining offsets. 

The questions of interest then are how important the low- 
frequency information contained in At is in reducing the 
the stripiness in the maps and how it depends on filter- 
ing out progressively higher frequencies. These questions 
can be rephrased as the following: suppose we observe 
two rings during time 2T. If there is no sky signal, how 
well can we constrain a relative offset, S, between the two 
rings using the data tl The variance in the offset obtained 
by generalized least squares is [i4W t _1 i] where now 
A = [0, 0, • • • , 0, 0, 1, 1, • • • , 1, 1]* describes the effects of the 
offset on the time stream. If the noise is stationary with 
power spectrum as in Eq. 1 then the (inverse) variance is 



^- 2 = iE 



I Wo (ATT 



(5) 



where the sum is over frequencies /k = k/2T from up to 
the Nyquist frequency and Wo/2 is the Fourier transform 
of the pointing matrix A which is given, up to a phase, by 



Wo (x) = j (irx) 



(6) 



where jo{x) = sin(x)/x is the spherical Bessel function 
of order zero. If the marginalization is performed the am- 
plitude of the inverse noise power spectrum is set to zero 
for all frequencies / < / mar g and as correspondingly in- 
creases. The rate of that increase depends on the assumed 
noise power spectrum parameters as shown in Fig. 1. It is 
clear that the mere presence of the I// component in the 
spectrum is nearly as harmful as the marginalization itself. 
That is because in such a case the noise power at frequen- 
cies / <C /knee is sufficiently high to effectively "weight 
out" (i.e., to marginalize over) the lowest frequency tail. 
Therefore in the absence of the I// noise the losses suf- 
fered due to the low-frequency marginalization are the 
largest. By the same token the time domain constraints are 
the strongest in a white noise case, though the dependence 
of a s on /knee disappears once / marg > / kn ee- For values 
of /marg of interest here (i.e., / marg < / spin < 1/2A S ) and 
an y /knee(> 0) the following approximate formula holds 
to within few percent, 



^r 2 ~ 2 

°s - °t 



,/kncc 



(7) 



log (/knee + /marg + /l) ~ log (/ ma rg + /l ) ' 

here, /i = I/2T ~ 3.5 • 10~ 5 Hz. This formula together 
with Fig. I are the main results of this Section. 

Note that in a more realistic case the available time 
stream data set will be composed of many ring segments of 
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Fig. 1. The dependence of the time domain constraint on rel- 
ative offsets as a function of the low frequency marginalization 
threshold (/marg) - Three curved lines correspond to the fiducial 
LFI, HFI and HFI-like white noise cases, respectively from top 
to bottom. For comparison, the nearly horizontal lines show 
the noise per beam-size pixel for LFI (top line) and HFI (bot- 
tom). The dotted vertical line shows the spin frequency of the 
planck satellite. 

length T. Therefore, rather than the single parameter case 
as in Eq. 5 we should consider a case with multiple offsets 
each assigned to a different stretch of the time stream and 
derive the constraints also for the relative offsets of the 
non-adjacent segments. Instead of deriving the appropri- 
ate formulas (which are a straightforward generalization 
of the single offset case considered above and follow the 
same basic steps as presented in the next Section), we just 
note that the relative offset of two segments becomes pro- 
gressively less constrained if the time separation of these 
two segments increases. That is a consequence of the ex- 
tra noise power present at the low frequency end of the 
spectrum either due to a generic, 1//, noise component 
or the low frequency marginalization. Therefore if there 
is any gain as a result of incorporating the time domain 
constraints on the relative offsets it should be the most 
significant on the shortest time scales. Hence we conclude 
that the formulas derived above, Eq. 5 & Eq. 7, indeed 
provide the best case answer to be confronted with the 
pixel domain constraints derived later on. 

5. Pixel (ring) domain considerations 

5.1. Ring- noise properties. 

Let A be now a pointing matrix from the time domain to 
a ring pixel domain. Then, as usual, noise correlations in 
the ring domain can be expressed as, 



K=(AW t - 1 A) \ 



(8) 



With crossing time of a single beam-size pixel denoted as 
irpix (= ispin/n-rpix), the pointing matrix, for a single ring, 
is given by: 



and is zero otherwise. Here we number the pixels in a way 
that ring pixel is observed at t = 0. Given a consecutive 
and linear pixel numbering along a ring circumference we 
have, 

A (t, k') =A(t+(k'- k) t ipix , k) , (10) 

for any two pixels k and k' . 

Rewriting Eq.(8) in a Fourier domain (and assuming 
circulancy of Aft) we have, 

-l 



(FA) 1 diag Pift)- 1 (FA)) 



here F is a Fourier operator. From Eq.(10) 
(FA) k = exp (-L2nf t k t rpix ) (FA) k=0 . 
From Eq.(9) 

(^)k=0 = W (ft irpix) Wl (ft ispin, T/t spin ) 

where 

sin (n + l)7ra; 



W x (x, n) 



■ exp (— i mrx), 



(11) 



(12) 



(13) 



(14) 



sin nx 

and Wo was defined in Eq. 6. Let us assume translational 
symmetry along the ring, neglecting therefore "switch- 
ing on and off" effects, the noise correlation in pixel do- 
main can be expressed via "ring domain" power spectrum, 
V (/ r ), where / r is a ring domain frequency expressed in 
Herz. The inverse noise spectrum in the ring domain is 
then given as, 



V(frY 



''spin 



= £ \Wi (ft t 

Afe, f t 

x P(/ t )- 1 exp(- t 2 7 r(/ t 



, T/t spin )\ 2 \W (f t t lpix )\ 2 



f r ) Afc f r pix) 



E 



\Wt (ft t apiD , T/t spin )\ 2 \W (f t t r 1 



^rpix/ 



(15) 



x W 1 x (n(f t - / r ) trpix, n rpix ) PC/t) 

Here, the superscript denotes a real part. The noise cor- 
relation matrix can be computed via a Fourier transform, 

Afr = F T diag [P (/,)] F. (16) 

Though this result looks very intuitive, the expression for 
V (/ r ) is less so. Nevertheless, for the set of parameter 
values as described in Sect. 3, the resulting ring-domain 
power spectrum turns out to be basically a time-domain 
spectrum resampled on a sparser grid of frequencies corre- 
sponding to the ring domain. The corresponding ring do- 
main /knee is nearly identical to that in the time domain, 
and the resampling ensures that the marginalization over 
all frequencies / < / ma rg does not cause any singularities 
of the ring domain noise correlation matrix in addition to 
that at zero frequency as long as / mar g < /spin, i-e., in that 
case V (/ r ) _1 > 0, if / r ^ 0, and V (f T ) is well-defined. For 
the beam-size pixels, we find that the off-diagonal noise 
correlations in the ring domain are usually less than a 
percent of the rms noise value per ring pixel, given by, 



A (t, k) = 1; t — t mod t rp i X = k t s 



(9) 



^rpix ^ 



(17) 
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In the following, wc will therefore neglect these off- 
diagonal correlations and consistently keep assuming that 
the number of independent, uncorrelatcd ring elements is 
equal to the number of beam size ring pixels, n rp j x . It is 
important to point out that the assumption made above 
does not automatically force the noise correlations of a fi- 
nal map to be also diagonal. We discuss that issue further 
below. 

Having estimated the noise correlation matrix, W r , we 
can also estimate the ring map using Eq. 4. In the follow- 
ing, in the spirit of the destriping methods we will recon- 
struct the sky maps via combining all the precomputcd 
sky rings paying a particular attention to the noise corre- 
lation patterns of the resulting sky maps and the precision 
of the recovery of the relative ring offsets. 

5.2. Offsets from ring crossings. 

In this Section we will proceed assuming that there are 
no off-diagonal noise correlations within each ring. This 
assumption has been justified above. We will show below 
that this assumption does not mean that there will be no 
correlations in the final map, made as a superposition of all 
the rings on the sky. We will estimate these correlations as 
well as dispersions and compare the results with the time 
domain constraints derived earlier. 

Yet again, let us define a simple map-making prob- 
lem, this time the projection is performed from sky pixel 
domain, as for example defined by healpix (Gorski et 
al. 1998), to the ring domain. From a set of all ring pixels 
for all the sky rings we choose only ring-pixels which cross 
with at least one other ring. We form a vector, r, of all es- 
timated temperature values for all such ring pixels. These 
are derived from the time ordered data as described in the 
previous Section. We assume that each of the ring pixels 
included above corresponds to a certain uniquely defined 
sky pixel. Though this may not always be the case, as 
some ring pixels may span more than a single sky pixel, 
it is a useful simplification which we expect to be broken 
only occasionally. We number the sky pixels consecutively 
counting each pixel once and form a vector of the corre- 
sponding true sky temperatures denoted, s. Consequently, 
any ring pixel temperature, n, can be modeled as a sky 
temperature in a corresponding sky pixel, sj , plus a ring- 
pixel instrument noise, rii and its own specific offset, x^. 
The latter contributes to any measurement made along a 
given ring. Hence we get, 



r = As + n + Bx = [A, B] 



+ n. 



(18) 



In the above equation, A is a pointing matrix assigning 
which sky pixel was observed at which measurement, and 
a matrix, B, is deciding to which ring a given ring-pixel 
belonged and therefore which offset x is to be added to 
it. Under our assumptions, the noise correlation matrix, 
W r = (nn*), is diagonal and all of its diagonal elements 
are identical and denoted as, ov 2 . 



Our task here is to solve this system of equations to 
estimate the noise properties of the recovered "map" , m, 
which includes estimates of both the actual sky, s, and 
ring offsets, x. The resulting noise correlation matrix is as 
usual given as, 



M s = {[A,B] t N- 1 [A,B}) 



A* A 
B l A 



A l B 
B l B 



(19) 



Hereafter we will denote the matrix on the right hand side 
of this equation with, A. This matrix has a block struc- 
ture with two diagonal blocks describing the correlations 
between the sky pixels and the offsets respectively, the two 
off-diagonal blocks characterize the cross-correlations be- 
tween them. Of course, not all ring pixels are involved in 
this map-making problem explicitly, but only those which 
are observed during at least two different ring scans. The 
correlations between remaining pixels are equal to the 
correlations between offsets of the rings they happen to 
be intersected by. For wcll-cross-linked scanning strate- 
gies, where most of the pixels are observed multiple times, 
from various directions, the pixel-pixel correlations need 
to be computed directly or based on considerations such 
as those presented here. In either case the level of stripi- 
ness in the maps can be measured as an RMS difference 
between the amplitudes of two rings given by, 



os = V 2 (ooff 2 - X oS ), 



(20) 



where cr ff and X g are the variance of the ring offsets 
recovery and a cross-correlation between two offsets. 

In general calculating Af s in Eq. 19 is difficult and com- 
monly requires extensive numerical computations, which 
easily become prohibitively expensive. Instead of doing 
that here we will consider a set of toy models to gain in- 
sight into the structure of the resulting noise correlations. 

5.2.1. two rings 

Let us begin with a pair of two intersecting rings (or more 
generally, basic scan paths). Let the number of crossings 
(i.e., pixels in common between the rings) be, n croS s. To 
break a degeneracy due to the unknown and unrecoverable 
absolute offset of these rings, we constrain the offset of 
the 1st ring to be 0. Recalling that every crossed pixel is 
observed twice, the matrix in Eq. 19 in this case reads, 



A = 



A 1 A, A l B 
B l A, B l B 



2 
2 



1 1 



1 n c 



(21) 



where the last column and row correspond to the offset of 
the second ring. The inverse of A is readily given by, 



A- 1 = 



A 1 A, A l B 
B l A, B l B 



1 



2n c 
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Across H~ 1 1 

1 Across 1 



^cross H~ 1 2 



• (22) 



The following observations can be now made: 

1. for a sky pixel, the variance decreases with the number 
of crossings very slowly and quickly reaches an asymp- 
totic value of cr r /-\/2; 

2. the off-diagonal terms for sky pixels are non-zero and 
arise as a result of solving for the unknown relative 
offsets between different rings. These correlations are 
however suppressed by a factor oc n cross + 1 with re- 
spect to the diagonal terms and therefore become pro- 
gressively less important as the number of crossings 
between two basic scan paths increases; 

3. the error in the determination of the relative offset 
between two rings is given by 

4. the correlations of the offset with the sky pixels in- 
creases accordingly so the off-diagonal term is always 
a half of the offset variance. 

For two special cases: one where two rings on the sky 
cross only in two points and the other when two rings fully 
overlap, we have, 

• two circles - two crossings (n cmss = 2): 



V3 



Coff = cr r ; 



(23) 



• two overlapping circles (n c 



^rpix 



cr x = er r , 



Coff = (T r 



2n T 



V2 ] 



< <7 X < cr r ; 



(24) 



X « 



Here, n rp i x is the number of independent ring-pixels along 
the circumference of each ring and in the case of PLANCK 
considered here it is given as n rp i x — 2000 ^> 1, justifying 
the inequalities and approximations made above. 

5.2.2. simple ring chains 

Let us focus on circular rings on the sky from now on. We 
will assume that those can either cross in two different 
ring-pixels or remain disjoint, neglecting therefore - with- 
out any loss of generality - cases of two tangent rings. Yet 
again let us begin with a pair of crossing rings. We can 
add a third ring in a number of ways. Clearly the most 



efficient way, from the view point of the offset constraints, 
is to add this ring in a way to make it cross each of the 
two former rings in two different pixels. The least efficient 
way, on the other hand, is when all three rings end up hav- 
ing a total of two pixels in common. We can keep adding 
more and more rings following these two prescriptions. In 
the following we analytically consider, in some detail, the 
resulting two extreme rings pattern on the sky. 
• «rin g rings all crossing each other in two fixed sky pix- 
els: 

A picture for this configuration can be that of scan paths 
following great circles on the sky intersecting only at the 
poles. Such a scanning strategy cannot be realized in prac- 
tice for a number of reasons (e.g., finite angular size of the 
detector focal plane, finite beam size) but it is a useful 
mathematical limit to consider. Again we force the offset 
of the first ring to be 0. The resulting correlations are then 
given by, 



A' 1 = 



A 1 A, 
B t A, 



A B 
B l B 



(25) 



^ring ~l~ 1 


faring 1 


1 


1 


1 


^Tiring 




~2 


~2 


~2 


^ring 1 


faring ~r~ 1 


1 


1 


1 


2?2 r j n g 


2?7. r i n g 


~2 


~2 


~2 


l 


l 


1 


1 


1 


~2 


~2 


2 


2 


1 

~2 


1 

~2 


1 

2 


1 : 





Here the left upper 2-by-2 block describes the correlations 
between the noise in the two crossing pixels, and the bot- 
tom right (n r i ng — l)-by-(n r i ng — 1) block shows the cor- 
relation and uncertainties in the recovery of the relative 
offsets in this case. Two things are evident: 

1. The noise per each of the two pixels very quickly 
reaches an asymptotic value of a r /\/2 when the num- 
ber of rings increases, and at the same time the 
two "antipodal" pixels become nearly 100% correlated 
(Janssen et al., 1996). 

2. The offsets variance and cross-correlations do not de- 
pend on the number of rings passing through the two 
pixels at all and they are equal to their respective val- 
ues as in the two-ring case considered before (Eq. 23). 

The stripes in a map produced for this scanning strategy 
would not go away with an increasing number of indepen- 
dent rings, as as = a T (Eq.20). 

• «rin g rings all crossing each other at two different sky 
pixels: 
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In this case a convenient way of visualizing a correspond- 
ing ring assembly is to think of a linear configuration with 
identical rings each shifted with respect to a previous 



"ring 

one by some small and constant displacement. The total 
number of all ring crossing is n cross = n ring (n r ; ng — 1) , and 
there — 1 relative offsets where, as before, we set 

the offset of the first ring to 0. The inverse of the matrix, 
A, can be then represented in a block form as follows, 



A- 1 = 



A' A, A t B 
B l A, B l B 



I 



2n r 



C D l E l 
D F G* 
E G H 



(26) 



where the meaning of the blocks is following: C describes 
the correlations between the pixels crossing involving ring 
0; F contains the correlations of all the other crossings and 
D defines the cross-correlations between these two sets of 
pixels. H describes the correlations between the offset of 
all rings, and E and G show their cross-correlations with 
corresponding sets of pixel crossings. Each of these ma- 
trices can be explicitly computed as it is shown in the 
Appendix. For the main purpose of this paper it is suffi- 
cient to notice that: 

1. In the limit of many rings (n r i ng £ 10) the variance 
of any crossed pixel is approximately a r /V2, and the 
cross-correlations are suppressed by a factor oc n r ; ng . 

2. The error in the reconstruction of the offsets is 
2cr r /^/n r i ng with their correlations equal to cr r /n r i ng 
(see Eqs. A. 2 & A. 5) and both decrease relative to the 
pixel noise variance, if the number of rings increases. 



Consequently, we have as = cr r \/2/n r ing. It is therefore 
apparent that this scan pattern not only better constrains 
the level of the map stripiness but also that it produces 
lower pixel-pixel cross-correlations than the scan consid- 
ered earlier. 



5.3. Getting (more) real. 

The considerations presented above are clearly idealized, 
though they capture a number of features of actual sky 
scans and therefore can serve as a guidance in understand- 
ing the properties and effects of scanning strategies. The 
major omission seems to be neglecting the non-zero width 
of the actual rings on the sky - a factor which becomes 
increasingly important for scans with a large number of 
rings. In such cases new effects can start playing a role 
due to additional overlaps between rings due to crowding. 
Therefore in the limit of large, n r ; ng , our conclusions may 
need amending. 

Let us consider the scanning strategy with two fixed 
crossing points with a large number of rings. Partial over- 
lap between different rings at high latitudes becomes un- 
avoidable leading to a quick improvement in the level of 
the map stripiness over the numbers quoted above. How 
rapid the improvement really is depends on the width of 
the rings and their numbers. For example if we consider 
rings of FWHM width and assume that there are no gaps 



left between the rings at the equator, we find that it cre- 
ates <~ 2 sin -1 (1/3) /7r n rp ; x ~ 0.2 n rp ; x of extra near 
crossings, i.e., ring pixel overlaps at which their centers 
are displaced relative to each other by less than FWHM/3. 
These crossings suppress the relative offsets by an extra 
factor oc ^/0.2 n rp j x ~ 20, making this scanning strategy 
more comparable to the other one from that point of view. 

For the scan with all rings crossing each other in two 
different pixels on the sky, there is a limitation on how 
many rings can be a part of such a configuration. If we 
assume again that all the rings have a width equal to the 
beam FWHM, an upper limit on the number of rings is 
given by the number of beam-size pixels which can be 
fitted along a quarter of the circumference of a ring, i.e., 
"ring ~ 5400 sin 6» /FWHM. In fact we find that for the 
PLANCK scanning pattern with an 85° opening angle, this 
limit can indeed be achieved within at most a factor of 2. 

It is important to note that as we assume that all rings 
are uncorrelated, it really does not matter if the rings used 
for producing a map were scanned by a single detector or 
many different detectors. As long as all the considered 
rings intersect each other in two different ring pixels, the 
upper limit on a number of rings remains unchanged, de- 
termining how well we can do in terms of constraining the 
relative ring offsets. As a result, combining the data of 
many such detectors can only improve the noise correla- 
tions of the coadded map inasmuch as the single detector 
maps were suboptimal. Moreover, once the limit of a maxi- 
mal number of intersecting rings has been reached, adding 
extra data will not bring any further gain either in terms 
of the resulting noise level per pixel or in terms of a lower 
relative level of off-diagonal correlations. 

By contrast, in the case of multiple detectors which 
perfectly follow the same trajectory on the sky, combin- 
ing the maps produced by each of them will lower the level 
of both diagonal and off-diagonal noise correlations by the 
same factor (oc n<j e t) if the noise properties of the detectors 
are identical). This results in lower noise in the final map 
but does not enhance the diagonality of the noise correla- 
tion matrix with respect to the maps made of the data of 
each of the detector separately. In particular such a pro- 
cedure will not suppress "stripes" , if those are present in 
the single detector maps. 



6. Comparison. 

6.1. Total intensity maps. 

In Sections 3 & 4 we have derived the limits on the offsets 
due to the time and pixel domain constraints respectively. 
Here we will compare them for the particular case of a 
PLANCK-like experiment. 

Let us begin with the time domain constraints. From 
Eq. 7 (see also Fig. 1) for the set of PLANCK-like values 
(Sect. 3) we find that erg ~ 20 - 70yuK (hfi and LFI re- 
spectively) if no marginalization is applied. That value 
becomes ~ 50 — 170^iK if the low temporal frequencies 
all the way up / sp i n ~ 16mHz arc marginalized over. It 



R. Stompor & M. White: Low frequency modes and map-making for PLANCK 



9 



is apparent that both these values are much too high to 
facilitate production of a high quality map of the sky. Yet 
the only robust way of improving on these numbers in the 
time domain is by decreasing the 1// part of the time do- 
main noise power spectrum. To get these numbers down to 
a micro-Kelvin level requires /knee of the order of a tenth 
of mili-Herz (Eq. 7), i.e., two orders of magnitude below 
the expectation for planck. 

The other way to alleviate the problem is to exploit 
the pixel domain constraints. That can be accomplished 
by a choice of scanning strategy. Given the results in 
Sections 5.2.2 & 5.3 we can estimate the level of stripiness 
for PLANCK-like scanning strategies, i.e., composed of re- 
peatedly scanned rings on the sky. For our fiducial value 
of a beam FWHM equal to 10' and any of the discussed 
scans, we get n r i ng ~ 500 and as ~ 1.7— 5/iK. These num- 
bers are well over an order of magnitude lower than the 
best (i.e., with no marginalization) time domain results, 
demonstrating that the stripes in PLANCK-like maps will 
be predominately constrained by the ring crossings. As 
discussed earlier a full optimal map- making procedure ap- 
plied to the time stream data divided into segments of the 
length much longer than T could utilize both types of the 
information, leading to an improved combined constraint 
on the ring relative offsets. However, given our estimates 
above, the resulting constraints can be tigthened only by 
at most 1% over the constraints derived from the pixel 
domain only. (We have allowed generously here for factors 
of order of few possibly missing in our pixel domain anal- 
ysis, and made use of a fact that both constraints should 
be combined in a noise weighted fashion.) 

That shows that the two-step map-making approach 
is a nearly lossless way of making the sky map from the 
PLANCK data. Moreover, again in the case of the optimal 
map-making, the marginalization over the long temporal 
modes, relaxing the time domain constraint by a factor 
~ 2.5, can lead only to an increase of ^ 2% of the average 
level of stripes in the final map. That demonstrates that 
for either choice of the map-making approach it is possible 
to produce the final sky map from PLANCK, nearly optimal 
and free of the parasitic contributions residing at the low 
temporal frequency end, ^, / sp m- 

The question, whether the achievable suppression of 
the variance of stripes down to a few micro-Kelvin level, 
as estimated here, is sufficient, will depend on a particular 
application the produced maps are to be used for and 
needs to be answered on a case-by-case basis. 

6.2. Polarized maps. 

So far we have focused only on the total intensity maps. 
However, some of the results can be generalized to give in- 
sight into polarized map-making with planck. This issue 
was discussed in some detail by Revenu et al.(2000), here 
we briefly highlight some of the related problems. 

Let us assume that we have four detectors following 
each other across the sky. Each detector is a total power 



detector but measures only photons with a specific po- 
larization as defined by a front mounted polarizer. For 
simplicity we will assume that all detectors sample the 
sky at precisely the same points along the scan trajec- 
tory and that the respective orientation of the polarizer 
is different in each of them, but fixed with respect to the 
scan direction. For definiteness in a numerical example 
below we assume these angles to be <pi = 0°, <p2 = 45°, 
03 = 90° and <j>4 — 135°, respectively. As total offsets of 
any of three Stokes parameters can not be recovered, we 
can set them to any arbitrary values and force the offsets 
of three data streams to be zero. The offset of the fourth 
detector stream is to be determined from the data. We de- 
rive the appropriate formulas in Appendix A. 2.1. We find 
that in the limit of many beam-size ring pixels the rms 
noise per pixel is twice higher for the Q and U maps than 
for the T map. Consequently the respective time domain 
constraints on the relative offsets of the two polarized ring 
maps has to be scaled accordingly, 







' i 






V2 






.V2_ 



where <j g is given by Eq. 5. The values in Eq. 27 are to 
be compared with the pixel domain constraint derived for 
each of the Stokes parameters. 

The derivation of the pixel domain constraint requires 
a conversion of the recovered Stokes parameters from the 
ring to global coordinates. For the latter we will use the 
usual coordinates where the preferred direction is set along 
the meridians on a sphere. In general such a transforma- 
tion introduces an extra angle dependence to the "pointing 
matrix" of the map-making problem including the offsets 
(Eq. 18). That is because what used to be just an offset 
in a ring coordinate system in the global coordinates may 
vary from a pixel to a pixel as a different rotation may 
be needed, to perform the transformation from the ring 
to the global coordinate system. As a result the crossings 
will constrain differently the relative offsets of two rings, 
depending on the precise geometry of the intersection. 

Although all of that makes the precise computation 
rather cumbersome, for the two examples of the scanning 
strategies considered here the general problem simplifies 
allowing for a fast estimate of the effect. In the case of 
scanning along great circles intersecting only at the poles 
the ring and global coordinate frames basically coincide 
and the derivation of the pixel domain constraints follows 
from our previous considerations leading to analogous re- 
sults (though with both types of the constraints rescaled 
as in Eq. 27). 

In the case of a scan pattern made of rings intersecting 
each other in two points, a conversion from the ring coor- 
dinates to global coordinates is necessary, requiring pre- 
cise knowledge of a particular scan geometry. However, it 
seems natural to expect that the cumulative loss with re- 
spect to the total intensity case should be ~ y/2 as a result 
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of averaging sine and cosine functions over different possi- 
ble angles between two rings at the crossing point. In fact 
this intuition can be made more quantitative and shown 
to be essentially correct (see Appendix A. 2. 2). Hence also 
in this case we derive the same results, within a factor of 
y/2, as in the total intensity case. 

Consequently our major conclusion about the irrele- 
vance of the low (sub-spin) frequencies for the sky map 
quality holds also in the polarization case. Note that again 
our statements are comparative and do not state that the 
residual level of stripiness is satisfactorily low, but rather 
that it is mostly determined by the pixel domain con- 
straints. 

7. Summary. 

We have investigated the importance of the low temporal 
frequency modes on a quality of the sky maps as antic- 
ipated from the planck satellite. We have shown that, 
for plausible scanning strategies and predicted instrumen- 
tal noise properties, marginalization over frequency modes 
lower than the satellite spin frequency does not cause any 
significant increase in the stripiness of the resulting opti- 
mal maps. This conclusion is based on a fact that the ma- 
jor constraint on the level of stripes is due to the overlap of 
the scan rings on the sky and not due to the low frequency 
modes contained in the time stream. Our results also sup- 
port the idea that the two-step map-making, where the 
maps of the scan rings are first computed from the time 
ordered data, and only then combined together to produce 
a map of the sky, is nearly optimal and therefore does not 
compromise the quality of the resulting map. Though we 
have neglected in our considerations of the scanning strat- 
egy possible complications, such as precession or nutation 
of the satellite spin axis, we expect that those are bound 
only to strengthen our conclusion as they commonly lead 
to an increase of a level of scan cross-linking. Therefore, 
in the context of planck the results derived above seem 
to be quite general and demonstrate the robustness of the 
mission design with respect to the long term systematic 
effects which may be present in the time ordered data. 

It is important to emphasize, that a major factor be- 
hind this result is the presence of 1// noise with an /knee 
frequency of the order of mili-Herz. In such a case the re- 
sulting high noise power at low frequencies leads to long 
term variations on its own, the variance of which exceeds 
the constraints imposed on the stripes by the scan cross- 
linking. If, however, the low frequency excess power were 
absent, what would require suppressing /knee by two or- 
ders of magnitude, both the time and pixel domain con- 
straints would be important in reducing the overall strip- 
iness of the resulting map (Mennella et al., 2002). 

We have also studied the properties of the noise cor- 
relations in pixel domain. We have shown that though 
for the beam-size pixels the noise properties within each 
ring on the sky are close to white noise, the final map 
composed of many rings can display a non-negligible level 
of off-diagonal pixel-pixel correlations. We have pointed 



out that though these can in principle be as high as a 
half of the pixel variance, for the two scanning strate- 
gies considered here the off-diagonal elements are typi- 
cally suppressed by a factor ~ 400 — 500, with respect 
to the pixel variance, and are therefore of a similar order 
as the pixel-pixel noise correlation within each ring due 
to the 1// noise. We therefore conclude that the pixel- 
pixel noise correlation matrix for PLANCK-like maps will 
be strongly diagonal dominated. Whether this observation 
justifies neglecting the off-diagonal terms in a statistical 
analysis of the planck maps may depend on a particular 
application. 
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Appendix A: Estimating pixel-pixel noise 
correlations. 

A.l. Total intensity case with every pair of rings 
crossing in two different pixels 

The symmetries in this case are somewhat less apparent 
than the case treated in the main text. First of all, we 
break the symmetry between rings on the sky by arbitrar- 
ily choosing the ring the offset of which is assumed to be 
zero. Also the difference between correlations of a pair of 
pixels which belong to the same ring and those which do 
not needs to be made adding yet another layer of struc- 
ture to the matrix. Nevertheless, a great deal of insight 
can still be gained without the need for massive (and in 
practice often prohibitive) numerical computations. 

Hereafter we elucidate block-by-block the structure of 
the matrix A as defined in Eq. 26. Let us start with the 
diagonal blocks. The matrix, C can be decomposed into 
2x2 block matrices, such that, 



C 



Cdiag Caff' 
Coff 



Coff 
Coff* 
^off 



c J 



c 



off 



c. 



(A.l) 



here the ith diagonal block express the correlations of ev- 
ery two pixels belonging to the ring and the ring i, while 
the off-diagonal blocks describe the cross-correlations be- 
tween these pairs of pixels. Note that all these pixels 
belong to the ring and hence their variance is some- 
what lower than all the other pixels considered below. The 
blocks can be written explicitly as follows, 



1 



1 



Cdiag — 



Coff = - 



i 

i i 
i i 



+ i 



(A.2) 



(A.3) 
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The matrix, F, encodes the correlation between all 
pairs of pixels belonging to the same (non-zero) ring with 
all other (non-zero) rings. Again, it can be represented as, 



-^diag ^off 



^off ' • -^diag . 

where the diagonal blocks read, 
3 3 

Tiring ~r 



(A.4) 



-^diag — 



•ring 

3 



(A.5) 



and the off-diagonal blocks come in two flavors depending 
on whether the two pairs of pixels in question belong to 
the same ring or not, 

5 

i r 



F oS = 



1 1 



(A.6) 



The correlation matrix for the offsets (H) is given as, 
'2 1 ••• 1" 



H 



1 



: 1 

'•• 1 
1 2 



(A.7) 



The off-diagonal blocks are somewhat more cumbersome 
as many more subcases need to be kept track of. To rep- 
resent cross-correlations between the rings offsets and the 
pixels belonging to the Oth ring, let us define first two 
vectors, 



E„ 



-1 



-1] 



(A. 



Now we can represent the matrix, E, as a vector of block 
matrices, each of the size (n r i ng — 1) x 2, i.e., 



E 



E, 



block 



block 



(A.9) 



where each of the blocks is made with (n r ; ng — 1)/2 rows of 
the 1st type above and (n r i ng + l)/2 rows of the 2nd type. 
The indices numbering the blocks in the above equation 
are just to distinguish different ordering of the basic vec- 
tors within a block. In a similar manner we can represent 
the G matrix, but this time the basic vectors need to be 
defined as, 



Gr, 



-1 



-l]x 



(A.10) 



All these calculations are clearly quite mundane and 
using these results directly in a case of particular ring ar- 
rangement on the sky may not be necessarily straightfor- 
ward. However, the results presented above demonstrate 
explicitly that one can gain some insight into the struc- 
ture of the correlations of PLANCK-like maps using sim- 
plified semi-analytical considerations. There are a number 



of important conclusions, which can be drawn from these 
results. Let us first recall that we have separated a fac- 
tor l/2/n r i n g in Eq. 26, which now needs to be taken into 
account. Thus we find that for this kind of strategy all off- 
diagonal terms scale as l/n r ; ng with a coefficient at most 
of the order of few. The same turns out to be true for the 
offsets, the variance of which, as expressed by the H ma- 
trix above, decreases with a number of rings in the scan 
in the same way as the level of the off-diagonal terms. 
However, the variances of the sky pixels themselves are 
nearly independent on the number of the rings for any 
large n r ; ng and approximately equal to a r /^/n r i ng . 



A. 2. Polarized case. 
A. 2.1. Single ring. 

Let us first consider a single ring on the sky observed by 
a four polarized detectors. As in Sect. 6.2 we assume that 
all the detector polarizers are at fixed angles with respect 
to the ring, which arc 0i = 0°, 02 = 45°, 03 = 90° and 
04 — 135° respectively. We also constrain the offsets of 
the three time streams to be zero and let the offset of the 
fourth detector (= 64) to be undetermined. As before we 
first compute a maximum likelihood ring-map for each of 
available time streams and then convert those into maps of 
the three Stokes parameters in the ring coordinate system. 
The former step is feasible as we assume that the orien- 
tations of the detector polarizers with respect to the ring 
direction does not change with time. Consequently all the 
results of Sect. 5.1 apply here directly. That latter step is 
again treated like a map-making procedure, however with 
the noise in the ring pixel domain assumed now to be diag- 
onal (Sect. 5.1). The corresponding pointing matrix needs 
to relate the measurements of any of the detectors with 
ring maps of the Stokes parameters. For a single ring pixel 
it would therefore read, 
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m 




r Q 


+ 






r\j 




ri3 




64 




Tl4 



(A.11) 



This equation is straightforwardly generalizable for the 
case of an arbitrary number of the ring pixels. For definite- 
ness hereafter we order the ring pixels either in a way that 
first n rp ; x values correspond to the total intensity ring- 
map, next n rp ; x values to the Q-map and [/-map for the 
Stokes parameters maps or simply concatenate together 
maps for each detector for the detector maps. The last 
(3n rp i x + 1) value is the relative offset of the fourth time 
stream. With all these assumptions, the noise correlation 
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matrix for the Stokes parameters determined in the ring 
coordinate frame reads, 



TV = 



Att Atq Atu Ats 
■Atq Aqq Aqu Aqs 
■Atu -Aqu Auu Atjs 
■Ats A vs A55 



*Q<5 ^US 

where the blocks are given by, 

n rp i x + 1 1 • • • 1 

1 : 



(A.12) 



A 



TT 



A 



TU 



: ••• •• 1 

1 • • • 1 n rpix + 1 

-2 -2 ••• -2" 



•• -2 

-2 -2 



(A.13) 



(A.14) 



-4qQ = 



2n r 



2n r 



(A.15) 





"2n rpix + 4 2 • 




2 








2 








A\j\j = 










(A.16) 










2 
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2 2n 


rpix "T" 4 ^ 






"-2" 




' -4" 






Ats = 




, A\js — 




, -4<5<5 = 


= [4]. (A.17) 




-2 




-4 







And all the remaining blocks are equal to zero. As already 
remarked in Sect. 6.2 the noise matrix is strongly diago- 
nally dominated in the limit of large n rp i x . Note also that 
the lack of symmetry between the Q and U Stokes param- 
eters, conspicuous in the above equations, is solely due to 
our choice of which detector offset is to be determined as 
a result of the procedure. In fact in the presented here 
case the fourth detector stream does not contribute to a 
constraint on the Q component, as can be seen by directly 
inverting the relation in Eq. A. 11 for our choice of the po- 
larimeters' directions (Revenu et al. 2000). Consequently, 
the unconstrained global offset of the fourth detector map, 



64, docs not cause any extra cross-correlations of the Q- 
map with itself or any other Stokes parameters. In fact 
the symmetry is restored in the limit of a large number of 
ring pixels as in this limit the cross-correlation due to that 
offset uncertainty (~ 2CT t / A /n rp ix) becomes irrelevant. 



A. 2. 2. Multiple rings. 

Let us assume henceforth that the number of ring-pixels 
is large and, consequently, the off-diagonal terms in 
Eqs. (A.12-A.17) are negligible. In this Section, we there- 
fore consider a set of uncorrelated rings of the three Stokes 
parameters each defined in a corresponding ring coordi- 
nate frame. The offsets of each ring maps are undeter- 
mined. As in Appendix A.l we assume that each pair of 
rings intersects in two independent ring pixels. Hereafter, 
we focus on the Q and U rings only as spatial rotations 
do not couple temperature to polarization. The discus- 
sion of the temperature case has been already given in 
Appendix A.l, below we show how one can gain some in- 
tuition about the results in the polarized case. We define 
a projection from the ring domain to the sky pixel domain 
as in Eq. 18. The ring pixel map is then expressed as, 



[A,B] 



+ n, 



(A.18) 



where, as before, only intersecting ring pixels are included. 
This time, however, both the ring pixel map, r, and the 
sky pixel map s, combine together maps of two Stokes 
parameters, Q and U. For future convenience, we define 
them as follows, 



Q12 
U w 
Q21 
U21 

Ql3 

Ui 3 



^3n r i n g 1-Hring 



0„ 

Un rl 



nngi u nng 



Q12 
U12 

Ql3 
Ul3 



Un r , 



-l,n ri „ g 
-l,n r ; ntI 



s? 



(A.19) 



where, for instance, Qij denotes the value of the Q Stokes 
parameter measured at the intersection of rings i and j 
during the scan of the ith ring and therefore expressed in 
the coordinate system of the ring i; Qij is then a value of 
Q component in the global coordinate system, and 6® is a 
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Q offset of the ith ring. In this convention, and (j, i) 
correspond to two different pixels where the ith and jth 
ring intersect. Note that we have n r ; ng relative offsets and 
(^cross + "cross) values of Q and U in sky pixel domain 
to be determined using 4n cross measurements. The noise, 
n, is assumed to be uncorrelatcd and its variance given 
by a T V2 (see Eq. A. 15 & A. 16). With these definitions at 
hand, the sub-matrix A is then block diagonal, 



A\ 
A\ 







(A.20) 



with 4-by-4 blocks given as, 



cos 2cnij sin 2ay 

— sin 2a;j cos 2a-^ 
cos2/3ij sin2/3ij 

— sin 2/3jj cos 2/3ij 



(A.21) 



Here ay (or /Jy) is an angle between the ring coordinate 
system of the ith ring and the global coordinate system 
at the sky pixel denoted (i, j) respectively). Hence, 

A\ performs two rotations from the global to the ith ring 
coordinate system at the two crossing points of the ith 
and jth rings. 

The elements of the block B are either or 1, and an 
element is 1 only if the particular crossing correspond- 
ing to the ith row of the A matrix is between a ring j and 
some other ring, and was measured during the scan of the 
jth ring. 

In principle, having determined the structure of both 
matrices, A and B, we can calculate the inverse (up to 
a factor) of the noise correlation matrix in sky pixel do- 
main, A, as defined in Eq. 19. However, in practice for 
a general scan pattern, considered here, the computation 
of A turns out to be quite mundane and requires more 
specific assumptions. Instead, as an example, we calculate 
only the diagonal block of A - one corresponding to the 
correlation matrix for the offsets, and therefore being of 
the most interest for the consideration contained in this 
paper. To do that first note that the corresponding block 
of the A matrix can be computed by partition and reads 
(Eq. 19), 



Ass 



B t B 



B t A(A t A) 1 A l B 



From Eq. A.21 we have, 



A* A = 27 n , 



crows s*. n crosK 



and also 

B l B — n r i ng / nr 



(A.22) 

(A.23) 
(A.24) 



from the definition of B in Eq. A. 18. Here, I n> , n stands for 
an n-by-n unit matrix. Also from Eq. A.21, we find that 
AA l is block-diagonal with the (i, j) 4-by-4 block given as, 



[AA* 



where, 



1 


Cii 




1 



c ji 

1 



(A.25) 



cos 2 (ay 
sin 2 (ay 



■Aj) 
Aj)- 



(A.26) 
(A.27) 

These formulas allow us to perform the matrix operations 
on the right hand side of Eq. A.22 rendering, 

-4 4ci2 

-4 



A' 1 

Ass 



n r 



4 Cln rlng 



4C12 



(A.28) 



4 C±n r m g ' ' "ring 1 

In the case of many rings, i.e., n r i ng S> 4 > |4c;j|, the ma- 
trix on the rhs of the last equation is nearly diagonal and 
we finally derive the approximate expression on the vari- 
ance of the relative ring offsets in the case of the polarized 
maps. This reads, 



as ~<T ti — <7T\ • (A.29) 



<Ting 
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